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ABSTRACT 



Context. CoRoT-2, the second planet-hosting star discovered by the CoRoT satellite, is a young and active star A total of 77 transits 
were observed for this system over a period of 135 days. 

Aims. Small modulations detected in the optical light curve of the planetary transits are used to study the position, size, intensity, and 
temporal evolution of the photospheric spots on the surface of the star that are occulted by the planetary disk. 

Methods. We apply a spot model to these variations and create a spot map of the stellar surface of CoRoT-2 within the transit band 
for every transit. From these maps, we estimate the stellar rotation period and obtain the longitudes of the spots in a reference frame 
rotating with the star Moreover, the spots temporal evolution is determined. This model achieves a spatial resolution of 2°. 
Results. Mapping of 392 spots vs. longitude indicates the presence of a region free of spots, close to the equator, reminiscent of the 
coronal holes observed on the Sun during periods of maximum activity. With this interpretation, the stellar rotation period within 
the transit latitudes of -14°6 ± 10° is obtained from the auto-correlation function of the time integrated spot flux deficit, yielding a 
rotation period of 4.48 days. With this period, the temporal evolution of the spot surface coverage in individual 20° longitude bins has 
periodicities ranging from 9 to 53 days with an average value of 31 ± 15 days. On the other hand, the longitude integrated spot flux, 
that is independent of the stellar rotation period, oscillates with a periodicity of 17.7 days, which false-alarm probability is ~ 3%. 
Conclusions. The rotation period of 4.48 days obtained here is shorter than the 4.54 days as derived from the out-of-transit light 
modulation. Since the transit data samples a region close to the stellar equator, while the period determined from out-of-transit data 
reflects the average rotation of the star, this is taken as an indication of a latitudinal differential rotation of about 3% or 0.042 rad/d. 

Key words, extra-solar planets; star spots; stellar magnetic activity 
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1. Introduction 

Rotational modulations of the light curves of stars have been in- 
terpreted as due to the presence of starspots on the stellar surface 
(e.g., Kron, 1947). Despite being first detected on red dwarfs, 
nowadays spot activity has been observed on all main-sequence 
late-type stars and is attributed to photospheric magnetic fields, 
by analogy with sunspots. By modelling the rotational modula- 
tion of the stellar flux, it is possible to infer the brightness and 
relative sizes of starspots. Moreover, long-term data series pro- 
vide information on stellar cycles, similar to the 1 1 - year solar cy- 
cle, spot longitude distri butions, and lifetimes (e.g., iBerdvuginal 



cie, spot longitude aistnbui 
I2005t IStrassmeieit I2009D . 

The spot modelling of the disk-integrated light modulation, 
however, only allows for a few large spots on the stellar surface. 
This changed with the discovery of the first planetary transits 
of HD 209458 (ICharbonneaui lloool: iHenrv et all I2OOOI) . When 
the planet eclipses the star, it may occult solar-like spots on the 
stellar surface, pro ducing detect able variations in the light curve 
during the transits. ISilval (I2003h developed a model to infer the 
properties of starspots from transit photometr y. This model was 
later applied to the light-curve of CoRoT-2 (ISilva-Valio et al.L 
l201ft) . the second planet-hosting star discovered by the CoRoT 
satellite. Here we give continuity to our analysis of the spots on 



* CoRoT is a space project operated by the French Space Agency, 
CNES, with participation of the Science Programme of ESA, 
ESTEC/RSSD, Austria, Belgium, Brazil, Germany, and Spain. 



CoRoT-2, following the spots in a reference frame that rotates 
with the star. 



This same star had its lo ng-term photometric out-of-transit 
data analysed by L anza et al.i (2009) . The authors concluded that 
the star had two preferred longitudes on opposite hemispheres, 
spot lifetimes of the order of a month, and derived a lower limit 
for the amplitude of the surface differential rotation of only 0.7% 
from the migration of the active longitudes. 

Furthermore, both the rotation modula tion and transits of 
the CoRoT-2 light-curve were modelled bv lHuberetaT](l2010h 
with alow angular resolution, reaching 15° during transits. Their 
modelling yields a rotation period within the transit band of 4.55 
days, while a longer period is suggested for the non-eclipsed 
latitudes giving a latitudinal differential rotation AO. ^ 0.1 
rad/d. The authors also find two preferred longitudes on opposite 
hemispheres. Here we focus on the modelling of the photomet- 
ric modulations produced by starspots during tran sits, aiming 
at a significantly higher longitude resolutio n than Huber et 'M 
(12010), thanks to the method introduced by ISilval (l2003h . This 
approach allows a better characterization of the properties of the 
individual starspots along the occulted band, detecting spots as 
small as ~ 0.2 planetary rad i i, i.e., 2° of longitude on the stel- 
lar surface (ISilva-Valio et al.l 1201 Oi see also Sect. |4]below) and 
using them as tracers of stellar rotation. 
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By tracking the m otion of sunspots across the disc of the 
Sun. ICar rington' ('1863') noticed that their sidereal angular veloc- 
ity could be approximated by; 



Transit 22 



Q = A - B cos^ 



(1) 



where A ^ U ^SS day"', B = 2°87 day"', and 6 is the co-latitude 
(iThomas & W eiss. 2008). According to this equation, the solar 
equator rotation period is 24.7 days, that is, faster than that of 
the poles of 30.8 days. Evidence of diffe rential rotation in ot her 
late-type stars has also been found (e.g., iBames etal I l200.5h . 

The differential rotation of distinct solar activity features 
is not the same, but depends on the features themselves. For 
instance, solar coronal holes, and their photospheric counter- 
parts represented by magnetic unipolar regions, usually exhibit 
a remarkably rigid rotation in comparison to active regions and 
sunspot groups. We find something reminiscent of this behaviour 
in the rotation of the activity features in the photosphere of 
CoRoT-2 because the active longitudes trace a pattern rotating 
much more rigidly than the individual spots (see Section|5]l. 

The next section provides an overview of CoRoT-2 observa- 
tions, whereas the following section presents the method used. 
Previous results are stated briefly in Section |4] Sections |5]and|6] 
describe the spots longitude and the spot temporal evolution, re- 
spectively, once the stellar rotation period is calculated. The next 
section explains how the stellar differential rotation is obtained, 
and finally Section[8]presents the discussion and conclusions. 



2. Observations of CoRoT-2 

The Convection Rotation and planetary Transits (CoRoT) 
space telescope, launched in December 2006, has discovered 
15 planets so far using the method of transits (see, e.g., 
http://exoplanet.eu). The second planet discovered by CoRoT 
orbits a young and active star similar to the Sun, now named 
CoRoT-2, of spectral type G7 V with mass and radius of 0.97 Mq 
and 0.902 Rq, respectively. 

CoRoT-2 was continuously monitored for about 135 days, 
during which 77 transits of its planet were observed with high 
temporal resolution (32 s). The first two transits were discarded 
because of the l ower tempora l resolu tion. The data were reduced 
as explained in lAlonso eta D dloosl) with a resulting rms of the 
points of the white light curve of 0.0006 in units of the relative 
flux of the out-of-transit data. The CoRoT white passband ranges 
fr om 300 to 1 100 nm an d its transmission profile can be found 
in lAuver^ne et al.l (12009 ^. together with a detailed description of 
the instrument, its operation, and the data processing pipeline. 

The basic assumptions and the parameters of the planet and 
its orbit used for the present modelling were taken from previous 
observations dAlonso et al l l2008h . and are hsted below: 

- Circular planetary orbit of radius 6.7 Rs, whe re Rg is the 
stellar radius, and period of 1.743 days; indeed. !Gillon et all 
12010) find that the orbit is slightly eccentric with e cos a> - 
-0.0029 + 0.0006, where e is the eccentricity and a> the argu- 
ment of periastron, but this small eccentricity does not sig- 
nificantly affect the results derived from a simple model with 
a circular orbit; 

- Orbital inclination angle of 87?84, implying a projection of 
the central transit chord at a stellar latitude of -14?627 (the 
Southern latitude is arbitrarily chosen); 

- The stellar equator and the planetary orbit are coplanar; 
this is supported by the measurement of the Rossiter- 
McLaughlin effect that yields a misalignment of the pro- 
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Fig. 1. Example of the fit to transit 22. Left: Star with modelled 
nine spots. Right: observed light curve of the transit (black) with 
the resulting light curve from the model (thick grey line). 



jected spin and orbital angular momenta of only 7° 
(Bouchy et al., 2008); 



+ 5° 



- Planet radius of Rj, - 0.172 R^; this value was estimated 
from the deepest transit, that is, the one with the minimum 
number of spots (ISilva-Valio et al.Ll2010l) . 



3. The Method 

The model used here simulates a star with quadratic limb dark- 
ening as a 2-D image such as the one shown in the left panel 
of Fi g ure [1] A detailed description of the model is given in 
Silva' (2003'). Applications of such a model are described in 
Silva-VaHo (2008) and Silva- Valio et al. (2010), the latter specif- 
ically for CoRoT-2. The model simulates the passage of the 
planet, a completely dark disk, in front of the spotted star The 
light curve is computed as the sum of the fluxes coming from all 
the pixels in the frame, for a sequence of planet positions along 
its orbit at intervals of 2 minutes. 

Unless otherwise stat e d, the planet and star parameters are 
taken from l Alonso et al.l (l2008h . Moreover, the model for the 
spots assumes that: 



- The spots are circular and described by three parameters: 
their radius (in units of planet radius, Rp), intensity (in units 
of stellar photosphere maximum intensity at disc centre, Ic), 
and longitude; 

- The spot latitude is the same as that of the central transit 
chord on the disk of the star, that is, -14°6; 

- The spot longitude is the stellar topocentric longitude, i.e., 
with zero longitude on the meridian at the centre of the 
disk of the star at mid-transit; an accuracy up to ~ 1° can 
be reached in the case of CoRoT-2 (cf., e.g., iWolter et al.L 
.2009.) : 

- The spots are constrained to longitudes +70°, to exclude the 
distortions on the ingress and egress branches of the transit 
profile that cannot be measured as accurately as those falling 
in the central part of t he transit because of th e steepness of 
the light variation (see lSilva-Valio et"an 1201 0); 

- The effect of foreshortening when a spot is close to the limb 
is taken into account. 

This model was applied to all transit light curves and the 
three parameters (radius, intensity, and longitude) for each spot 
were sought simultaneously by minimizing the between the 
model light curve and the data. All fits were performed using the 
AMOEBA routine (Press et al., 1992). Initial guesses of the pa- 
rameters were obtain ed through the genetic algorithm PIKAIA 
(ICharbonneaulll995h . 



A. Silva-Valio and A. F. Lanza: Evolution and rotation of starspots on CoRoT-2 



Number of spots 



40 



30 T 



20 r 



Fig. 2. Histogram of the minimum number of spots per transit as 
derived by the method described in Sect. |3] 



The minimum number of spots per transit, M, needed to fit 
the Hght curve within the data uncertainty, was estimated as fol- 
lows. A model with M spots has 3M free parameters, because 
each spot corresponds to 3 free parameters (longitude, radius, 
and intensity). Therefore, the number of degrees of freedom 
of the model is s - N - 3M, where is the number of data 
points per transit (usually N - 211). Now, the problem is to find 
the minimum number of spots that yields an adequate best fit 
for each transit according to the criterion for hypothesis testing 
given by, e.g., Lampton et al. ( 1976). We compute the minimum 
value of the for each model with a different M starting from 
the case with a single spot (M = 1). A model with M spots is re- 
jected at a significance level y if its minimum chi square x^iM) 
exceeds x^iy), i-C-, the y-point of the;^'^ distribution defined by: 



Xoo 
J(y) 



(2) 



where / is the density of the probability distribution of the 
X^ with s degrees of freedom. By increasing M, the minimum 
X^(M) decreases until for a certain value of M, say M, we have 
X^{M) < x^^iy)- The minimum number of^£9ts_giving a n ac cept- 
able model is therefore M. Following iLampton et alT ( ri976l) . we 
assume a significance level y = 0.1 and compute thex^ cumula- 
tive distribution function in Eq. ^ by means of the IDL function 
CHISQR_PDF. 

Figure |2] shows the distribution of the best number of spots 
per transit, that is, the minimum number of spots that makes the 
fit statistically acceptable. A total of 392 spots were modelled 
during the whole period of observations (77 transits). The aver- 
age number of spots per transit was 5 + 1 . 

In principle, our modelling approach can suffer from the ef- 
fects of possible systematic errors, but we believe that their im- 
pact is negligible because: a) the model has enough degrees of 
freedom to account for a highly complex flux modulation dur- 
ing transits by increasing the number of small spots; b) the effect 
of bright faculae is not detected in our transit profiles because 
we never observe dips below the reference profile as derived 
from the deepest transit; c) systematic errors arising from limb- 
darkening parameters or the adopted relative radius of the planet 
are minimized by considering only the central part of each tran- 
sit profile (+70°). In any case, our choice of a confidence level 
7 = 0.1 excludes too small x^ values that may imply an over- 
fitting of the data or enhance the impact of any residual system- 
atic error 
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Fig. 3. Histograms of the spot parameters obtained from the fits 
to the 77 transits: spot intensity (in units of Ic) (left), spot radius 
(in units of Rp) (middle), and flux deficit due to spot (right). 



4. Spots physical parameters 

The location in time of a "bump" along the transit light curve 
provides unambiguous information on the longitude of the cor- 
responding spot. When the spot latitude is fixed a priori, the ra- 
dius and the intensity of the spot may show some degree of de- 
generacy, but in the case of a highly accurate light curve, such as 
for CoRoT-2, they are independently constrained by the duration 
and the height of the corres ponding burnp, resp ectively. This is 
confirmed by the analysis ofl Wolter et al.l (120091) who considered 
also spots centred at latitudes different from that of the central 
transit chord, thus studying the effect of this additional degree of 
freedom in the model. 

A quantity proportional to the blockage of stellar emission 
due to the spot, called relative flux deficit, can be defined as 
D - {I - fi) , where fi and S are the spot intensity (in units of 
Ic) and radius (in units of Rp), respectively. The deficit is propor- 
tional to 1 - / because darker spots are the predominant ones. In 
fact, a value of /, = 1 means no spot at all, since there is no con- 
trast with respect to the photosphere. The advantage of using the 
flux deficit, instead of the values of f, and S , is that it is virtually 
independent of the degeneracy between these parameters. 

Figure |3] shows the results of the spot intensity, radius, and 
flux deficit for the spots detected during all transits, a total 
of 392 spots. These are different from the results reported in 
Silva-Valio et al. (2010), where a maximum of 9 spots was as- 
sumed. Here, the number of spots fitted was smaller, about half 
than previously used, owing to the criterion applied to find the 
optimal number of spots (cf. Sect.[3]l. The new model gives mean 
values of the spot intensity (Figure |3] left panel) and the radius 
(middle panel) of (0.45 + 0.25) and (0.53 + 0.18) Rp, respec- 
tively. The corresponding distribution of the relative flux deficit 
is shown in the right panel of Figure |3]and has an average value 
of 0. 14+0.06. Because fewer spots per transit were considered in 
th e new fits, these average values are larger than those reported 
in lSilva-Vahoetan dloiol) . 

Simulations for a single central spot with different radius and 
intensity show that only the spots with a flux deficit greater than 
approximately 0.02 produce a signal in the light curve above the 
data noise level of 0.0006. Therefore, only the spots with flux 
deficit above 0.02 are considered in the following analysis and 
have their best-fit parameters plotted in Figure |3] Note that this 
selection criterion implies a threshold of 0.20 Rp for the spot 
radius and a maximum value of 0.92 Ic for the spot intensity. 

The longitudes of the spots obtained from the best fit to each 
transit are shown in Figure H where the size of the symbols rep- 
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Fig. 4. Top: Diagram of the spots longitude, size, and intensity 
for each transit. Darker spots are pictured in darker shades of 
grey. Bottom: Spots flux deficit integrated in time as a function 
of spot longitude. 

resents the relative diameter of the spots, and their colour their 
intensity. Black spots are those with the highest contrast. Only 
spots within -70° and +70° longitude and flux deficit greater 
than 0.02 are plotted. The bottom panel shows a histogram of 
the flux deficit vs. longitudes for these spots. 

5. Spot longitudes and stellar rotation 

The longitude of a spot obtained from the fit is the topocentric 
longitude on the star measured in an inertial reference frame. 
To convert this value into a longitude measured in a reference 
frame that rotates with the star, that is, fixed to the stellar surface , 
one needs the rotation period of the star. lAlonso et al. I (l2008h 
obtained a value of 4 .54 days for the mean stellar rotation period, 
whereas Lanza et al.1 (12009) found 4.522 + 0.024 days from the 
motion of the active longitudes. 

The rotational longitude, ySrot, is determined from the 
topocentric longitude, yStopo, as: 

A-ot=Aopo-360°^^, (3) 

"star 

where foA = L743 day is the planet orbital period, /"star the 
stellar rotation period, and n the transit number The rotational 
longitude is limited to the range +180° and coincides with the 
topocentric longitude at the mid-time of the first transit, i.e., 
when « = at HJD 2454242.7666. 

5. 1 . Comparison with tlie out-of-transit spot modeliing 

Before determining the best rotation period within the tran- 
sit band, it is interesting to compare the longitudes of the 
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Fig. 5. Spots longitude versus time for spots with flux deficit 
larger than 0. 1 are shown as white circles overplotted on the map 
of the spot filling factor obtained by Lanza etal. (2009) from 
out-of-transit data (their Fig. 4). Different colours indicate dif- 
ferent filling factors, with yellow indicating the maximum spot 
occupancy and dark blue the minimum. The bridges connecting 
active longitudes are labelled with A, B, and C (see the text). 



spots ob tained here with the out-of-transit spot maps as de- 
rived by iLanza etal] ( |2009) . First, the spot topocentric longi- 
tudes shown in Figure|4]were converted to rotational longitudes 
for a stellar period of 4.522 days. Fixing the first transit at HJD 
2454242.7666, the shift of the origin of the longitudes between 
our reference frame and that of iLanza et al.l (12009) is -1 16°03. 
Thus, our longitudes were shifted by -116 ?03 to overp l ot bot h 
results. Figure |5] is the same as Fig. 4 of iLanza et al. I (l2009l) . 
with the white circles representing the spots from the transit fits 
with a flux deficit greater than O.L The radius of each circle is 
proportional to the flux deficit of the corresponding spot. The 
agreement between the results of the two different modelling ap- 
proaches is remarkable. 

Some of the small and medium spots seem to trace small 
bridges between the active longitudes, even if they are not per- 
fectly superimposed with the Maximum Entropy bridges marked 
with A, B, and C in Figure. |5] If interpreted as the effect of mov- 
ing spots, these bridges can be used t o derive the differential ro- 
tation (cf., e.g., Frohlich et al., 2009). However, the transit maps 
suggest that they are not single migrating spots but instead are lo- 
cations where spots appear and disappear between the two active 
longitudes, like "hot spots" where activity is preferentially lo- 
cated between the active longitudes. Of course, not all the spots 
found from transit mapping are reproduced in the out-of-transit 
maps because short-lived spots do not produce a significant ro- 
tational modulation signal, while the main features found from 
the out-of-transit data are also found by transit mapping, includ- 
ing some evidence of those small bridges. Moreover, the spatial 
resolution of the transit mapping is higher than that of the out- 
of-transit modelling which is limited to approximately 30° -50°. 
Therefore, the spots mapped during transits show a more discrete 
distribution vs. longitude. 

5.2. Rotation period at tlie transit iatitude 

The rotation period of the transit band needs not be equal to the 
star average rotation period, because the star may present differ- 
ential rotation. To determine the rotation period of the stellar sur- 
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Fig. 6. Spot rotational longitudes for 6 different rotation periods 
(reported at the top of each panel) for the spots of all transits 
with flux deficit larger than 0.1. The size of the circles is pro- 
portional to the spot radius whereas the colour refers to the spots 
flux deficit. The light grey bands in the background represent 
the stellar longitudes visible to the observer Bottom: Spots flux 
deficit integrated over time as a function of spot longitude. 



face around the transit band, which covers -14°6 + 10° latitude, 
spot maps of the stellar surface within the transit latitudes were 
constructed, similar to those of Figure |4] Periods ranging from 
4 to 6 days were tried, with 0.01 day intervals. Figure |6] shows 
such maps for the periods of 4.30, 4.45, 4.48, 4.52, 4.54, and 
4.60 days. These maps show the relative size and flux deficits of 
the spots as a function of their rotational longitude (horizontal 
axis) and transit time (vertical axis). 

At first glance, one sees that the spots follow diagonal stripes 
in these panels, giving the impression that they are wrapped 
around the star too tightly or too loosely. That is to say, the stellar 
periods used are not appropriate for the latitude of the belt oc- 
culted during transits. A visual inspection of these panels shows 
that it is as if the spots of the initial diagonal stripes (slanted to 
the left) are being straightened out to an almost vertical tempo- 
ral evolution (for Ps,ar around 4.54 days) and then the diagonal 
inclination changes sign (slanted to the right) for periods longer 
than 4.54 days. Similar to adjusting the image of an old TV set. 

If the stellar rotation period at -14°6 latitude is taken as 4.54 
days, this would imply that CoRoT-2 has almost no differential 
rotation since this practical ly coincides with t he average period 
of the star as measured bv iLanza et al.l (l2009l) . The same result 
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Fig. 7. Left: Auto-correlation function of the spot flux deficit vs. 
longitude for four different stellar rotation periods as indicated 
by the different curves. Right: Width of the main peak at half 
the maximum value (FWHM) of the auto-correlation function 
versus stellar rotation period. 

has been obtained bv lHuber et al.l (1201 0^ as they derived the ro- 
tation period by minimizing the drift of the active longitudes. 

Nevertheless, a different interpretation is possible by analogy 
with the rotation of solar coronal holes. Specifically, the blank 
range between 0° and 40° longitude, free of spots (see Figure|6), 
can be interpreted as analogous with a solar coronal hole that 
extends to -14'?6 latitude and lasts for over 120 days. From solar 
studies, it is known that coronal holes rotate more rigidly than 
sunspot groups, and with an equatorial rotation period similar to 
the average solar rotation period (cf. Sect. |8). If we consider this 
to be the case here, we can trace the rotation of individual spots 
as they move across the active longitudes and the blank range. 

Therefore, concentrating on individual spots, by visual in- 
spection it may be possible to identify several lines of beads. 
Because of the large quantity of spots, we used only those with 
a flux deficit greater than 0.1 for this purpose. The possibility 
of tracing the longitude of individual spots vs. time is limited 
by the fact that a homogeneous longitudinal coverage of the 
occulted band requires at least five transits (cf., e.g.. Fig. 2 of 
iHuber et al.. 2009), i.e., ^ 9 days during which spots gener- 
ally evolve remarkably both in radius and intensity. Therefore, 
a quantitative criterion to find the best rotation period is to look 
at the auto-correlation function of the time integrated flux deficit 
as a function of longitude, for the most promising periods. They 
correspond to the maps with vertical strings of spots, i.e., show- 
ing the minimum drift of the individual spots vs. time. 

Figure |7] (left panel) shows the auto-correlation functions for 
periods of 4.45, 4.48, 4.54 and 4.55 days. The right panel of 
Figure Q displays the full width at half maximum (FWHM) of 
the main peak of the auto-correlation function versus the rotation 
period from 4.4 to 4.6 days. As can be seen from the figure, the 
period that yields the narrowest auto-correlation function is that 
corresponding to 4.48 d, indicated by the vertical dashed line. 

Therefore, hereafter the adopted stellar rotation period in the 
latitude band occulted by the planet is 4.48 days. A map of the 
surface of the star for each transit as a function of the rotational 
longitude for such a rotation period is plotted in Figure |8](upper 
panel). 

6. Spot temporal evolution 

After the stellar rotation period and thereafter the rotational lon- 
gitudes of the spots have been determined, it is possible to esti- 
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Fig. 8. Surface map (top) for spots with flux deficit larger than 
0.02 and time integrated flux deficit of the spots versus longitude 
(bottom) for a stellar rotation period of 4.48 days. The light grey 
bands in the background of the top panels represent the stellar 
longitudes visible to the observer 



mate the spot temporal evolution. Here the rotational longitudes 
are calculated by considering the star rotation period within the 
occulted latitude band to be 4.48 days. First, the surface of the 
star is divided into 20° longitude bins. Then the total flux deficit 
within each bin is calculated as a function of transit time. The 
resulting flux deficit in each bin is plotted in Figure |9] where a 
running mean of every five data points has been used to warrant 
a homogeneous sampling of all the longitudes. 

The temporal evolution of the spot surface coverage within 
each longitude bin is estimated by computing the Lomb-Scargle 
periodogram of the time profiles of Figure |9]with a maximum 
false-alarm probability of 0.01 (99% significance). The periods 
obtained from the time series analysis of the data range between 
9 and 53 days, with a mean value of 31 + 15 days. 

The time profile of the spot flux deficit integrated over all 
rotational longitudes is plotted in the top panel of Figure [TOl and 
is independent of the adopted stellar rotation period. A running 
mean every 5 points is overploted as a thick line. A periodic- 
ity is clearly seen in this integrated variation. Its Lomb-Scargle 
periodogram is shown in the bottom panel of Figure [10] and a 
prominent peak is seen at 17.7 ± 2.3 days. Coincidentally, this 
period is close to 10 times the orbital period of the planet, i.e., 
17.43 days, or approximately 4 times the stellar rotation period 
in the occulted band, i.e, 17.92 days. 

To determine whether the 17.7 d periodicity is real, we per- 
formed a random permutation of the transits, and recalculated 
the periodogram using the IDL SCARGLE routine with a maxi- 
mum period of 100 days. A total of 10,000 permutations were 
done and the histogram of the resulting periods is plotted in 
Figure [TT] As can be seen from this plot, the probability of a 
chance occurrence of the periodicity at 17.7 days (marked by 
the dashed line) is ~ 3%, therefore, we conclude that this peri- 
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Fig. 10. Top: Longitude-integrated flux deficit of all the spots in 
each transit. The thick line represents a running mean every 5 
transits. Bottom: Power spectrum of the total flux deficit plotted 
in the top panel. 



odicity should not be an artefact of the observational sampling, 
even though it is very close to 10 times the orbital period. 

The periodogram of the variation of the total spotted area as 
derived from the modelling of the out-of-transit light curve by 
Lanza et al. (2009) is shown in their Figure 7. It shows a highly 
significant periodicity at 29 + 4 days, whereas the third peak in 
order of decreasing power corresponds to a period of 8.7 days, 
i.e., to the first harmonic of the periodicity of the total spotted 
area in the occulted band. Although its power is significantly 
below the level corresponding to a false-alarm probability of 1 %, 
a component of the spot area modulation wit h that period canno t 
be excluded. The independent modelling bv lHuber et al.l (1201 Oh 
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Fig. 11. Histogram of the periods calculated from the longitude- 
integrated total flux deficit when the order of transits is shuffled. 
The dashed line at P=17.33 days represents the periodicity of the 
observed time series. 



Fig. 12. Differential rotation curve versus latitude for CoRoT-2 
(solid line) and the Sun (dotted line). The diamond shows the 
value of the rotation period of the transit latitudes. 



confirms the variat ion of the total spotted area vs. time found by 
iLanza etaTI (l2009l) . giving support to their results. 



7. Differential rotation 

The analysis of the out-of- transit data y i elds a stellar rotation pe- 
riod of 4.54 days (lAlonso et al.L l2008h . iLanza et all (l2009l) ob- 
tain a slightly shorter period of 4.522 days. While the out-of- 
transit period is an average over the entire latitude range cov- 
ered by the starspots, the period from our spot mapping refers 
to a restricted band of the stellar surface, that is, only the lati- 
tudes probed by the planetary transit, -14°.6 + 10°. Therefore, 
the period of 4.48 d found here refers to a latitude close to the 
equator and implies a differential rotation. This differential rota- 
tion could not have been detected if the dominant region free of 
spots, analogous with a solar coronal hole, had not been ignored. 

The Sun rotates differentially with an angular velocity law 
given by Eq. ([TJ. This angular velocity can be converted into a 
period as a function of co-latitude, 6: 



Psun = 



360° 



A-B cos2 e 



days, 



(4) 



where A and B are the same constants as in Eq. ([T). To obtain 
the equivalent law for the differential rotation on CoRoT-2, a 
function similar to Eq. © is used. This assumption is su pported 
by the Doppler imaging observations discussed by Barnes et alj 
(l2005h that show that a solar-like differential rotation law is gen- 
erally appropriate to fit the variation of the angular velocity vs. 
latitude in late-type stars. In the case of CoRoT-2, the values of A 
and B are determined by matching the period of 4.48 days where 
the co-latitude is - 90° - a - 75.4°, and 4.54 days where this 
function is averaged over all values of 6. These two conditions 
yield A = (80°5 + 0°3) day"' and B = (2?9 + 0°8) day-', thus: 



360° 



!oRoT-2 



80.5 -2.9cos2(90° - a) 



days, 



(5) 



where a is the latitude. The differential rotation of CoRoT-2 
given by Eq. (|5]) is plotted in Figure [T2] as a function of the 
latitude. For comparison, the differential rotation of the Sun is 
plotted in the figure as a dashed line (scale on the right side of 
the plot). The diamond marks the data point obtained from the 
transit spot mapping, that is, the CoRoT-2 period of 4.48 days at 
a = -14°6. 



8. Discussion and conclusions 

A total of 77 transits of the planet around CoRoT-2 were anal- 
ysed and the variations detected on the transit light curves were 
modelled as being due to flux deficits caused by circular spots 
having all the same latitude of -14°6, on the stellar surface. The 
physical cha racteristics of these spo ts were presented in a pre- 
vious paper ( Silva-Valio et al.L 120 10*). The spot parameters were 
recalculated here, yielding fewer spots, however the qualitative 
spots characteristics have not changed. Here we concentrate on 
the temporal evolution of the spots. 

Before the spot temporal behaviour can be estimated, it 
is necessary to know the rotatio n period of the st ar Out-of- 
transit data y ield periods of 4.54 (l Alonso et all |2008) or 4.522 
(iLanza et all 12009) days. The agreement of the spots longitude 
and flux deficit obtained from out-of-transit and transit mapping 
is remarkable (see Figure|5), and provides an independent confir- 
mation of the out-of-transit maximum entropy mapping method. 

Here, the goal was to obtain the rotation period within the 
transit latitudes from the spot longitude map. At first sight, it ap- 
peared that the surface rotation period of the transit latitudes was 
equal to the out-of-transit value, i.e., 4.54 days, implying that the 
star displays no differential rotation. However, if one interprets 
the lack of spots, or active regions, within the longitude range 
between 0° and 40° as the photospheric counterpart of a coronal 
hole extending to the equatorial regions (in analogy with the so- 
lar case), then it is no surprise that this region presents a rotation 
period similar to the mean value of 4.54 days. 

In the Sun, extended unipolar magnetic regions are of- 
ten found in the photosphere below coronal holes and share 
their rotational properties with them. Specifically, the prop- 
erties of solar coronal holes have been inv estigated by. 



e.g.. IShelke & Pandd (IT985[) . Ilnslev et al.1 (Il995h . iBumba et al.1 
(Il995h . I Ikhsanov & Ivano 3 (119991). There are basically two 
classes of coronal holes, those that are confined to the polar 
regions and dominate during periods of minimum activity, and 
those that appear at mid-latitudes or even in the equatorial re- 
gion, and predominate during solar maximum. The p olar coronal 
holes may also present equator-ward extensions (Insle v et al.L 
[1995). Insleyetal. (1995) studied the diff'erential rotation of 
coronal holes, and found that they rotate more rigidly than the 
photosphere. The mid-latitude coronal holes have an equatorial 
rotation period of 27.1 days, whereas the rotation of the equator 
ward extensions of polar coronal holes is even more rigid, show- 
ing almost no differential rotation. The rotation period at the 
equator of this class of coronal holes is 28.3 days (Ilnslev et al.L 
Il995i) . The equatorial rotation period of the mid-latitude coro- 
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nal holes is very close to the average rotation period of the solar 
photosphere. 

Ignoring the longitude range free of spots in CoRoT-2, we 
searched for vertical strings of starspots in the maps of Figure|6] 
while also inspecting the auto-correlation function of the spot 
time integrated flux deficit vs. longitude. According to this cri- 
terion, the best rotation period was found to be 4.48 days that 
minimized the drift in longitude of the individual spots and was 
taken as the rotation period of the stellar surface within the 
-14°6 + 10° latitude band. Using this rotation period for the star, 
the stellar differential rotation profile, shown in Figure [12] was 
defined. This profile shows that close to the equator the star ro- 
tates faster, with a period of 4.47 + 0.01 d, while the longest 
period of 4.61 +0.03 days is found at the poles, as obtained from 
Eq.dSL 

iLanza et al ] (I2009h found a lower limit to the amplitude of 
the relative differential rotation of ~ 0.7 - 1 % for this star, sig- 
nificantly smaller than the Sun. The relative differential rotation 
of the Sun is 22%. Here we obtain (f poie - Peq)/P = 3%, where 
P is the mean rotation period. In CoRoT-2 the out-of-transit light 
modulation is dominated by the lack of spots within the coronal 
hole close to the equator, yielding almost no differ ential rotation. 

Fro m a spot model with three spots only, iFrohlich et al.l 
(l2009l) find a mean differential rotation of 0. 1 1 rad/d and a max- 
imum differential rotation of 0.129 rad/d. Using the values ob- 
tained here, we find AQ = Qeq - flpoie - 0.042 rad/d . This is an 
intermediate v a lue be tween that found bv lLanza et aL l(l2009h and 
iFrohhch et al ] (I2009h . Considering the results of iBarnes et alJ 
( 20051) . ttie differential rotation of CoRoT-2 is approximately 
half the pole-equator angular velocity difference generally found 
in late-type stars of similar effective temperature and rotation pe- 
riod, suggesting that its spots are concentrated in an equatorial 
belt, as in the Sun. 

The temporal evolution of the spot surface coverage within 
20° longitude bins can be estimated. The average value of 
the spot flu x temporal variation is estimated as 31 ± 15 days. 
iLanza et"al] ([2009), from out-of-transit data, find that the total 
spotted area of CoRoT-2 oscillates with a period of 29 ± 4 days, 
close to the mean spot flux deficit periodicity estimated here. 

Another interesting feature is the temporal evolution of the 
average flux deficit, i.e., integrated over longitude in the occulted 
band. It shows a periodicity of 17.7 + 2.3 days with a false-alarm 
probability of ~ 3 %. This is very close to 10 times the orbital 
period of the planet or approximately 4 times the stellar rotation 
period found within the transit band. 

These results may suggest some kind of magnetic interaction 
between the star and the planet. Circumstantial evidence for star- 
planet interaction in the magnetic activity of the photosphere has 
been reviewed by Lanza ( 2008, 2009) and Shkolnik et al. (2009) . 
A major problem is the lack of a common and well-defined phe- 
nomenology in all the stars suspected to show such an interac- 
tion. Therefore, our suggestion should be regarded only as a ten- 
tative interpretation of the periodicities found for the total spot- 
ted area. 

An alternative interpretation may be that proposed by [Lou! 
(I2OO O') to account for the short-term cycles observed in the Sun 
close to some of the maxi_ma of the 11 -yr cycles, t he so-called 
Rieger cycles (iRieger et al.L Il984l) . For instance, lOliver et"an 
(Il998h detected oscillations with a period of ~ 158 days i n the 
total s unspot area close to the maximum of solar cycle 19. lLoul 
(I2OO O') interprets such oscillations as signatures of hydromag- 
netic Rossby-type waves trapped in the subphotospheric layers 
of the star. In the case of the Sun, the main periodicity corre- 
sponds to ~ 6 equatorial rotation periods. Since the rotation pe- 



riod of CoRoT-2 is only 4.5 days, this would correspond to a 
wave period of ~ 27 days, close to the periodicity of 28.9 days 
found by Lanza et al. (2009) from the modelling of the out-of- 
transit light modulation. Moreover, other periodicities, ranging 
from 2 to 6 rotation periods, have been re ported from the analy- 
sis of different solar data sets (lLoull2000l) . In the case of CoRoT- 
2, we found in Sect. [6] a periodicity of 17.7 days in the total 
spotted area of the occulted band, very close to 4 stellar rotation 
periods at that latitude, which would ten tativel y correspond to 
the 102-day periodicity reported by Loul (l2000l) for the Sun. Its 
first harmonic, corresponding to 2 rotation periods, could also be 
present in the Sun as well as in the variation of the total spotted 
area of CoRoT-2 derived from its out-of-transit light modulation 
(cf. Sect. [6]). 
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